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ABSTRACT 


This paper reviews defects in the diamond-type lattice with special 
emphasis on silicon which is used as the basic material in the fabri- 
cation of solar cells in space programs. The defects are classified as 
first, second, and higher generation defects, according to their com- 
plexity. Some of these defects include vacancies and interstitials (as 
in lithium, for example) which are first generation defects capable of 
migrating above certain temperatures and trapping charge carriers. 
These first generation defects, because of their mobility, can form 
more complicated defects, such as the Si-A and Si-E centers. The 
properties of such defects have been systematically tabulated. A short 
discussion is given on some simple models of the kinetics of such de- 
fects. A comprehensive bibliography is given at the end of the review 
which should help the uninitiated commencing research in this field. 
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STRUCTURE AND KINETICS OF DEFECTS IN SILICON 


by 

Martin M. Sokoloski 
Goddard Space Flight Center 


REVIEW OF DEFECTS IN THE DIAMOND-TYPE LATTICE 
Review of the Diamond Lattice 

The space lattice of diamond consists of two interpenetrating face-centered sublattices cen- 
tered at (0, 0, 0) and (1/4, 1/4, 1/4), respectively, as shown in Figure 1. Since the diamond structure 
is very inefficient in packing atoms as compared with the face-centered and body-centered lattices, 
it is a relatively empty structure. Each atom 
in the diamond structure is at the center of a 
tetrahedron with the nearest neighbor atoms at 
the vertices. Because of this, the inter-atomic 
forces are not only central but angular- 
dependent. This complicates matters when cal- 
culating the lattice vibrational spectrum (Refer- 
ence 2). A modified OPW method employing a 
pseudopotential has been successful in calcu- 
lating the band structure of elements crystal- 
lizing in the diamond structure (References 3 
through 6). Carbon, silicon, and germanium are 
some of these elements, the latter two being 
elemental semiconductors. Currently, silicon 
is being extensively used as solar cell material 
in the space program. Therefore, most atten- 
tion will be focused on it. 

First-Generation Defects in the Diamond Lattice 

First-generation defects are defined here as simple defects which are composed of not more 
than one defect. 



Figure 1 -Crystal structure of diamond, showing 
the tetrahedral bond arrangement (after W. 
Shockley, Reference 1). 
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The Lattice Vacancy 

One of the simplest types of defect is the lattice vacancy— an atom missing from its lattice site 
shown as A in Figure 2. The electronic rebonding after the defect has relaxed is shown in Fig- 
ure 3. There are three equivalent bond positions for the model; (26, 34), (23, 64) and (24, 63). A 
study in the temperature range of 14° to 20° K indicates an activation energy for reorientation of 
bonds from one equivalent configuration to another of only 0.01 to 0.02 eV at 20° K. This reorien- 
tation can also be accomplished by applying a uniaxial compressional stress along a cubic axis 
(Reference 7). The vacancy in silicon is capable of trapping electrons and existing in several 
charge states. In p-type Si, it is usually neutral, and in n-type Si a single or double negative state 
may exist depending upon the position of the Fermi level (see Figure 4). A positive-charge state 
exists but can only be seen if the sample is illuminated, since the trapping level of the positive 
vacancy usually lies below the acceptor impurity level. The vacancy is also capable of migrating 
throughout the lattice at room temperatures and below (see Figure 5). This is not surprising 




W n\W 

>VrV) n\y 

cfvWVY^o 

rf V/ Vf V V/ V/ V> 

VVAA) c /W 

crWWW\) 


A: LATTICE VACANCY 
B: HOST INTERSTITIAL 
C: IMPURITY INTERSTITIAL 


D: "INTERSTITIAL" OXYGEN ATOM 
E: DONOR IMPURITY IN SUBSTITUTIONAL POSITION 
F: ACCEPTOR IMPURITY IN SUBSTITUTIONAL POSITION 
G: DIVACANCY 


Figure 2— Schematic illustration of some defects. 
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since, as has been noted, the lattice is relatively open. Currently, it appears that the migration 
energy of the vacancy is dependent on its charge state, being 0.33 eV for a neutral vacancy and 
less than 0.20 eV for either negative-charge state. At this point it still is not known if the 0.20 eV 
energy is associated with a single or double negative state (Reference 10). 

The Interstitial 

Another simple type of defect is the interstitial shown as B and C in Figure 2. B represents 
the host atom in a nonlattice position. In the case of silicon, this would be the silicon atom. Thus 
far, the interstitial silicon atom has escaped detection (Reference 7). All experimental evidence 
indicates that the interstitial silicon is very mobile— much more than the vacancy. C represents 
an impurity interstitial which also can migrate throughout the lattice. Perhaps the most investi- 



O ATOMS ABOVE PLANE 
OF PAPER BY 0.072 a 0 



(a) TETRAHEDRAL SITE - NEAREST NEIGHBOR AT 0.433 a 0 ; 
NEXT NEAREST NEIGHBOR AT 0.500 a 0 . 


(b) HEXAGONAL SITE - NEAREST NEIGHBORS AT 0.415 a 0 ; 
NEXT NEAREST NEIGHBORS AT 0.649 a Q . 

THE PLANE OF THE PAPER IS NORMAL TO THE 
111 DIRECTION. 

Figure 6- Interstitial sites in the diamond lattice (after 
K. Weiser, Reference 11). 


gated impurity interstitials have been gold, sil- 
ver, and lithium, of which lithium is of most 
interest. The interstitial can occupy two differ- 
ent symmetry sites in the lattice, the interstitial 
position halfway between two next neighboring 
atoms along a (001) lattice axis, called the tetra- 
hedral configuration, and the position halfway 
between two next neighboring tetrahedral inter- 
stitial positions called the hexagonal configura- 
tion (see Figure 6). Most interstitials carry an 
electric charge when they enter the lattice. 
Lithium, copper, silver, and gold do so in silicon 
(References 12 through 16). However, the last 
three impurities can also exist as substitutional 
defects, and their diffusion is by the "dissocia- 
tive mechanism" (Reference 17). These elements 
can dissolve in Si or Ge in two states, inter sti- 
tially and substitutionally. To diffuse substitu- 
tionally, the host crystal must have a supply of 
vacancies available. Therefore, this facet of 
diffusion will depend strongly on the dislocation 
density in the crystal. The experimental values 
for the migration energy of Li in Si and Ge are 
0.66 and 0.51 eV, respectively (Reference 18). 

It is also known in radiation damage studies 
that interstitial Si atoms can replace dopant 
atoms, the dopant atoms becoming interstitials. 
For instance, an EPR spectrum arising from A1 
has been observed in aluminum-doped Si (p-type) 
irradiated at 4° K. There is evidence that other 
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group m impurities behave similarly (Refer- 
ence 8). The trapping levels of the interstitial 
A1 atom are shown in Figure 4. 

The Si-O-Si Configuration 

This defect, which is usually seen in crys- 
tals grown in quartz crucibles, is shown as D in 
Figure 2. The defect will hereinafter be called 
the tT interstitial tf oxygen atom. This basic de- 
fect was first observed by its 9jLt-band excita- 
tion about v 3 shown in Figure 7. An 11.7/I- 
bandis associated with the corresponding defect 
in Ge (Reference 19). There are six equivalent 
positions; however, Hrostowski and Alder point 
out that the barrier for reorientation from one 
equivalent position to another is so small that at room temperature the Si-O-Si molecule is essen- 
tially a free rotator (Reference 20). If a uniaxial stress is applied at an elevated temperature, the 
<111> direction will become nonequivalent and the oxygen atom will jump to a completely new 
orientation. The energy required to do this is 2.561 a 0.005 eV. This is equivalent to the energy 
needed for the diffusion of oxygen in the silicon lattice (References 21 and 22). The energy re- 
quired for migration of the defect is so large that treating the oxygen atom as a stable defect at 
room temperature and below is a good approximation. This defect plays an important role in 
second-generation defects, discussed below. 

Substitutional Defects 

Substitutional defects are shown as E and F in Figure 2. These defects consist of foreign 
atoms occupying a lattice site of the host crystal. These impurities may enter as neutral or 
charged substitutional defects, E being a group V donor and F being a group m acceptor. Although, 
when in a substitutional position, these atoms look like Si, nevertheless they are defects and should 
be treated as such. This will be more evident in the following section. 

Second-Generation Defects 

Early in the history of defects in Si and Ge, all trapping levels were thought to be due to iso- 
lated vacancies and interstitials, first- gene ration defects. This model soon proved to be incom- 
plete. In reality, the vacancy not only can be a trapping and recombination center, but can migrate 
throughout the lattice and interact with other defects to form second-generation defects. 

The Divacancy 

In this defect configuration, bonding occurs between a-d and a-d’ with the unpaired electron (in 
the solid line) between b-b f , as shown in Figure 8. Besides this, there are two other equivalent 
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configurations: bonding a-b and a 9 -b f with the 
unpaired electron on d-d T , and bonding b-d and 
b^d* with the unpaired electron on a-a f (see 
Figure 9). By stressing the crystal uniaxially, 
the majority of defects can be preferentially 
aligned in one direction. These preferentially 
aligned defects can then be frozen in at a low 
temperature. Then by a slow isochronal anneal, 
the energy of reorientation of the electronic 
bond can be measured. In the temperature range 
from 40° to 110° K, line broadening is observed 
in the EPR spectrum of the defect (Reference 24). 
The energy of activation for electronic reorder- 
ing is 0.06 eV. Nevertheless, another degree of 
freedom of the divacancy is about a line joining 
of this axis, an energy of 1.3 eV is required 
(Reference 24). This is actually the energy needed for the migration of such a defect. These vari- 
ous energies are classified in Figure 10. The trapping levels associated with the divacancy are 
rather complicated, as can be seen from Figure 4. The energy needed for the complete breakup 
of the divacancy is approximately 1.85 eV; i.e., the sum of the energy needed to separate the di- 
vacancy to the first nearest neighbor, 1.3 eV, and the corresponding energy needed to move the 
vacancies away from each other more than one nearest-neighbor distance, 0.65 eV (Reference 24). 
During low-energy electron irradiation, most of the divacancies are correlated; this is demon- 
strated by the anisotropy, R , of the vacancy orientation versus bombarding electron energy, where 
R is the ratio of the number of divacancies orientated along the beam axis to that along another, 
<111>, direction. For bombardment energies less than 1.0 MeV, R = 2.0 (Reference 25). At these 
bombardment energies, the rate of growth of divacancies is less governed by the diffusion of va- 
cancies than at higher energies. The 1.8ju, 3.3j u and 3.9/j, absorption bands have since been asso- 
ciated with electronic energy transitions of the divacancy (Reference 26). These bands are 


DooJ 


tin] 



Figure 8— Model of the divacancy, 
the vacancies. In order to change the orientation 



Figure 9—Three equivalent electronic orientations of the divacancy 
(after R.P.I. progress report. Reference 23). 
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associated with excitation charge carriers into 
traps. As a result there is no photoconductivity 
connected with these absorption bands. 

The Oxygen- Vacancy Defect in Si or the 
Si-A Center 

In migrating throughout the lattice, the va- 
cancy may encounter the structure in Figure 7. 
If so, it becomes trapped by the "interstitial" 
oxygen atom and forms the oxygen-vacancy 
complex shown in Figure 11. The configuration 
has an acceptor level at E c - 0.17 eV. The 
oxygen-vacancy defect can capture an electron. 
After this, a hole can be captured and recom- 
bined with the electron. The defect is neutral 
when capturing an electron, and negatively 
charged before it captures a hole. The capture 
cross section is for electrons, cr e ~ 10” 16 cm 2 , 
and for holes, a h ~ 10” 14 cm 2 (Reference 27). 
Again, there are six equivalent configurations 
for the oxygen-vacancy defect, the energy re- 
quired for electronic reorientation being 
0.20 ± 0.03 eV (Reference 28). Also, the defect 
itself can change positions corresponding to an 
applied stress. The energy required corre- 
sponds to 0.38 ± 0.04 eV for a neutral oxygen- 
vacancy defect (Reference 28). This energy is 
the thermal-activated switching of the oxygen 
atom from one pair of silicon atoms to another. 
The oxygen- vacancy defect gives rise to an im- 
purity vibrational band at 12 p (Reference 29). 
The trapping level at E c - 0.17 eV and the 12 ji 
band disappear upon an isochronal anneal to 
500° C, indicating an activation energy of 1.3 eV 
for possible breakup (Reference 30). 



i i i i r p_] i 

1 ST 2 nd 3 rd 4 th 1 V™ n + 1 TH 

VACANCY -VACANCY SEPARATION (ATOM SITES) 


NEIGHBOR POSITIONS ARE SHOWN IN ORDER OF 
INCREASING SEPARATION; THESE ARE NOT NECESSARILY 
IN A STRAIGHT LINE. 

Figure 10— The potential energy (in eV) of one neutral 
vacancy in the vicinity of another (after G. D. Watkins 
and J. W. Corbett, Reference 24). 



Figure 11— A model of the oxygen-vacancy complex or 
the Si-A center. 


The Phosphorus -Vacancy Defect in Si or the Si-E Center 

As in the oxygen -vacancy defect, the vacancy, if encountering a phosphorus atom in a sub- 
stitutional position, can become trapped and thus form the phosphorus-vacancy defect shown in Fig- 
ure 12. Again there are three equivalent electronic orientations. The energy needed to activate 
the hopping of an unpaired electron from one configuration to another is 0.06 eV (Reference 31). 
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The energy required to change the direction of the phosphorus-vacancy axis by rearranging atoms 
through a compression along a 110 axis at high temperature is 0.93 ± 0.05 eV (Reference 31). This 
reorientation amounts to diffusion of the phosphorus- vacancy pair through the lattice. It is not 
known if the phosphorus -vacancy defect breaks up or diffuses as an entity to traps, as pointed out 
by Watkins. Some annealing experiments give 0.93 eV as the activation energy for annealing, 
suggesting diffusion to traps (Reference 27). An argument by Watkins and Corbett gives the energy 
of annealing as 1.18 eV. This information is condensed in Figure 13. This defect has a deep ac- 
ceptor level at E c - 0.40 eV. So far, there is no evidence of an associated lattice vibrational band. 
The defect has a hole capture cross section of 1.1 x 10^ 13 cm 2 (References 27 and 32). 



Figure 12— A model of the phosphorus-vacancy complex 
or the Si-E center. 



PHOS. -VACANCY SEPARATION ( ATOM SITES) 


NEIGHBOR POSITIONS ARE SHOWN IN ORDER OF 
INCREASING SEPARATION; THESE ARE NOT NECESSARILY 
IN A STRAIGHT LINE. 

Figure 13— The potential energy (in eV) of a vacancy in 
the vicinity of a substitutional phosphorus atom (after 
J. W. Corbett, G. D. Watkins, R. M. Chrenko, and 
R. S. McDonald, Reference 29). 


Other Impurity -Vacancy Defects 

The silicon crystal can also be doped with other group V elements besides phosphorus. Studies 
on silicon doped with As, Sb, and Bi indicate that each of these impurities can form an impurity- 
vacancy defect in Si with trapping levels around E c - 0.40 eV and configurations similar to the 
phosphorus-vacancy complex. An empirical formula has been obtained for the activation energy of 
annealing these centers: 


E a = 3. 6r - 3.0 eV , 

where E a is the activation energy for annealing the defect, r is the tetrahedral covalent radius of 
the impurity atom of group V in A° (Reference 33). 

The A s - vacancy defect has a hole capture cross section of 4 x 10" 14 cm 2 (Reference 27). 

There is also a possibility of group DI atoms contributing to impurity- vacancy defects. For 
instance, in A1 -doped vacuum floating- zone Si, an EPR spectrum is observed to grow; this is 
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associated with the disappearance of the free vacancy spectrum along with the A1 (Reference 27) 
hyperfine interaction (Reference 8). No trapping levels have so far been associated with this Al- 
vacancy defect. 

Ion-Pairing 

Perhaps the most extensively studied interstitial in the Si lattice is the Li + ion which enters 
the lattice as a donor (References 12 and 34). In p-type Si, the Li + interstitials are known to inter- 
act with the p-type dopant atoms to form relatively immobile Li-dopant ion pairs (References 35, 

36, and 37). It is also known that the Li + interstitial will interact with the "interstitial" oxygen 
atom to form (Li-0) + defects. The dissociation or annealing energy of such complexes is obtained 
from the ratio of the effective drift mobility of Li in samples with oxygen concentrations up to 
1.3 x 10 18 cm" 3 to known drift mobility of free Li, and is 0.42 eV (Reference 37). The absorption 
peaks for Li-B pairing have been observed by W. G. Spitzer and M. Waldner and are reproduced 
in Table 1 (Reference 38). Some of their conclusions are: 

1. The modes at 564 or 584 and at 655 or 681 cm" 1 do not involve much Li motion. 

2. The position of these modes does depend upon a change of the local environment around the 
boron atom due to the Li ion, paired with the boron. 

3. The force constant between the paired Li and B is not great enough to cause a characteristic 
Li-B vibration. 


Table 1 

Number of Dominant Frequency Absorption Bands Upon "Annealing” 
the Oxygen- Vacancy Defect in Li. 


Isotopic pair/cm 1 

1J B &7 Li 

10 B & 1 i 

n B & 6 Li 

10 D & 6 
B L i 

Low-frequency peak 

522 

522 

534 

534 

Central-frequency peak 

564 

584 

564 

584 

High-frequency peak 

655 

681 

657 

683 

Peak due to 





impaired boron 

620 

644 




Higher-Generation Defects 

It was mentioned that the phosphor us- vacancy defect could migrate at an energy and become 
trapped, possibly forming another defect. Also, the migration energy of the oxygen- vacancy com- 
plex is such that this entity could migrate and form higher-order defects. Obviously the classifi- 
cation of these defects along with other more complicated defects can get out of hand. Therefore, 
the classification of defect conglomerations composed of three or more first-generation defects 
will be called higher -gene ration defects. 
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CWM Model for Annealing the Oxygen- Vacancy Defect 

Upon monitoring the decay of the impurity vibrational band of the oxygen- vacancy defect in Si 
at 829 cm” 1 (12p.), Corbett, Watkins, and McDonald noticed the appearance of other new bands at 
887 cm” 1 (11.215/i), 904 cm" 1 (11.004/x), 968 cm -1 (10.276/x), and 1000 cm-i (9.948 /i). The rise 

of the 887 cm' 1 line is interpreted as the dif- 
fusion of an oxygen-vacancy defect to an "inter- 
stitial tT oxygen site as shown in (a) of Figure 14, 
and is hereinafter called the CWM I defect. The 
remaining three bands would be caused by the 
addition of another "interstitial" oxygen atom to 
CWM I. This configuration shown as (b) in Fig- 
ure 14 will be called the CWM 2 defect (Refer- 
ence 39). As stated in Reference 39, these 
models are highly speculative. It is not known 
if either CWM 1 or 2 acts as trapping or re- 
combination centers. Since the experiment of 
CWM, other spectra shown in Table 2 have been 
seen upon annealing the oxygen- vacancy defect. 
Four of the lines have been seen and correlated 
by two or more investigators, including the 12/i, 
11.3/1, ll.Oju, and 9.0 /i lines (References 40, 41 and 42). (The slight differences between the values 
of Corbett, et al. and those of the other three investigators are due to a different measuring tem- 
perature.) This clearly demonstrates that the oxygen- vacancy defect need not decompose into its 
first-generation defects, but may form higher-generation defects. This may also be the case with 
other second-generation defects upon annealing. 



Figure 14— Tentative models for (a) the 887 cm -1 ir band 
and (b) 904-, 968-, and 1000-cm" 1 bands (after E. M. 
Pell, Reference 35). 


Table 2 

Associated Spectra Due to Li and Boron Pairing. 





Dominant Band in (Wave Numbers, 

Cm' 1 ) in Microns 



Investigator 

Cm -1 

(488) 

(834) 

(865) 

(894) 

(904) 

(922) 

(932) 

(936) 

(968) 

(984) 

(1000) 

(1106) 


20.5 

12 

11.6 

11.3 

11 

10.8 

10.74 

10.69 

10.28 

10.17 

9.95 

9 

R. E. Whan 
(Reference 40) 


X 

X 

X 


X 

X 

X 




X 

J. C. Corelli, et al. 
(Reference 41) 



X 

X 









X 

J. W. Corbett, et al. 
(Reference 39) 



X 


X 

X 




X 

X 

X 

X 

H. Y. Fan 
(Reference 42) 

X 

X 










X 


10 



The Silicon K Center 


M 


This defect is one of the predominant para- 
magnetic defects introduced into p-type silicon. 

This center requires oxygen and acts as a hole 
trap at E v + 0.30 eV. Its charge state can also 
change, depending on the position of the Fermi 
level (Reference 44). A speculative model pro- 
posed from EPR measurements by a group at 
RCA is shown in Figure 15 (Reference 45). Two 
oxygen-vacancy complexes can randomly mi- 
grate throughout the crystal and, upon a chance 
meeting, form such a defect configuration. The 
probability of this event occurring is more 
likely than two vacancies being trapped by a Si- 
O-O-Si molecule. Unfortunately, no impurity 

lattice vibrational mode has yet been seen or measured corresponding to this center. It would be 
interesting to see if any of the absorption peaks of the K-center correspond to any of those for 
CWM 1 or 2. 


Figure 15— RCA tentative model for the K-center. 



Conclusions 

All the kinetic properties of the various defects are given in Table 3. The contributors to the 
various experimental values which are unreferenced in the table can be found in the main text, and 


Table 3 

Defects in Li and Their Properties. 



Energy Parameters in eV 

Associated 


Defect 

Electronic 

Reorientation 

Migration 

Energy 

Annealing 

Energy 

Lattice 
Absorption 
Band (ju) 

Impurity Electronic 
Energy States 

Vacancy 

0.01 - 0.02 

0.30 (neutral) 
.20 (charged ) 0 

- 

- 

three charged states 
see Figure 4 

Silicon interstitial 

- 

0.02 

- 

- 

? 

’Interstitial" oxygen 

0.02 

2.56 

- 

9 M 

none 






double neg. 

Divacancy 

0.06 

1.3 

1.85 b 

? 

single neg. see 

neutral Figure 4 






single pos. 

Oxygen- vacancy 

0.20 0.03 

0.3F 

1.3 

12 n 

single neg at E c - 0.17 eV 

Phosphorus- vacancy 

0.06 d 

0.93 0.05 

1.18 e 

1.38 f 

? 

single neg at E c - 0.40 eV 

Li + 

- 

0.655 

- 

- 

acts as a donor 

(Li - 0) + 

- 

- 

0.42 

? 

acts as a donor 

Li + + B" 

- 

Continued: 

- 

see text 

- 

CWM 1 

- 

- 

- 

11.215 (887 cm" 1 ) 

? 

CWM 2 




9.948 (1000 cm-*) 

7 




10.276 (968 cm- 1 ) 
11.004 (904 cm' 1 ) 



K-center 

- 

- 

1 

i 

E v +0.3 eV, (hole trap) 


(a) May be a single or double negative charge state. 

(b) This is presumably for a neutral di vacancy. This energy probably depends on the charge state of the defect. 

(c) Energy required to break an Si-O bond for neutral Si -A defect 

(d) The energy required for the hopping of one unpaired electron between the three equivalent directions. 

(e) From EPR measurements. 

(f) From electrical measurements. 
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need not be repeated in the table. This list of defects is by no means exhaustive. Future develop- 
ments will add to the list and disagreements will persist, i.e. 

1. Quenching experiments in p-type Si report an energy of migration of the interstitial of 
0.31 eV (Reference 46). 

2. No theoretical calculations substantiate any values for the experimental values of the energy 
of migration of either the vacancy or interstitial (References 47 and 40). 

The vacancy is the defect which is responsible for most of the damage in the silicon lattice. With- 
out it, most second- and higher-generation defects would not exist. Therefore, the importance of 
establishing a well based theoretical calculation which agrees with the experimental value of the 
migration energy of the vacancy cannot be overstated. 


KINETICS OF DEFECT PRODUCTION AND MIGRATION 

The simple laws of mass action can be used with some success in following the various inter- 
actions of defects in solids. In this simple notation, the creation of an impurity-vacancy defect by 
a source of radiation, 0 , can be schematically written in the form of three reactions: 

4 > — V V + i , (1) 

K 

V + i annihilation , (2) 

V + I C , (3) 

where 

v is the free-vacancy concentration, 
i is the free-interstitial concentration, 

I is the free-impurity concentration, 

C is the 2nd-generation defect concentration. 

These concentrations are measured in atomic fractions. The side conditions are: 

i = V - C , 

I = l 0 - C , 

where i 0 is the impurity concentration at time t = t 0 . The reaction rates, k., are defined as 
follows: 


K i " v i exp(-E./kT) , 
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where 


is the coordination number, 
v. is the frequency factor, 

E. is the energy associated with the particular reaction, 
k is Boltzmann’s constant, 

T is the temperature. 

This model assumes, of course, that the interstitials are captured at sinks with rates much less 
than the K. f s defined. Actually, the reaction rate for trapping of interstitials by sinks such as the 
surface is so small compared to the other reaction rates that it can be neglected. Reaction 1 is 
just the creation of first-order defects (Frenkel pairs) through some source term <£, with an 
efficiency factor K 0 . Reaction 2 is the annihilation of these first-generation defects by chance 
meeting. Finally, reaction 3 is the trapping of vacancies by impurities to form second-generation 
defects, in this case impurity-vacancy defects. The differential equations become: 

dV 9 / X ( 4 ) 

"dF = -IW + K 3 C - KjV 2 + (k 2 - Kj) VC + K 0 4> 


dt" = iW ' K 3 C " K 2 VC • W 

These have been solved by computer program; asymptotic analytic solutions have been obtained for 
long times, low temperatures, and K 3 ^0. 

Reaction 3 is shown as being reversible. However, in the discussion of higher -generation 
defects, it is shown that the production of such defects is not only possible but likely. In annealing 
studies where the source term is now zero, the kinetic scheme would have to be replaced by more 
complicated reactions. For example: 


K, 

C 3 1 

V + I 

(6) 



c + c 

4* C 

(7) 


* 2 * 

C + V — 

K. 

-^ C 3 • 

(8) 

and, if any interstitials are still present, 



K, 

V + i * 

ann i h i 1 a t* i on : 


otherwise K. = 0. The differential equations would then be 


u 

CO 

1 

II 

CJ|^ 

- k 4 c 2 - k s cv , 

0) 
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( 10 ) 


dc 2 


dt ~ 

k 4 c 2 

dC 3 

dt “ 

k 5 cv 


dV 

dF = K 3 C - K,V. 


(ID 

( 12 ) 


In the case of the oxygen-vacancy defect in Si, reaction 7 would possibly be the production of the 
RCA K-center. At the present, higher-generation defects associated with the impurity oxygen 
atom have only been seen, but these reactions should not be ruled out for other impurities, such 
as the phosphorus-vacancy defect. In the case of measuring some electrical parameter, the 
"annealing" of the sample is usually seen as the trapping or expulsion of trapped carriers. If, 
upon "annealing", the resultant higher-generation defects were capable of trapping additional 


carriers, then there would be "reverse annealing 1 

TOTAL DOSE SAMPLE 



Figure 16— Reverse annealing of electrical parameters 
in Si (after T. Nakano, K. Nakasima, and Y. Inuishi, 
Reference 48). 


r . This can best be seen by counting charges: 

V + I - K 2 ► C - Me* 

C + — „C 2 -- (M + ( ) e* . 

The result of the anneal will be the additional 
loss of 1 charges. However, an isochronal 
measurement of the intensity of the defect ab- 
sorption bands together with photoconductivity 
will only show forward annealing of c and pro- 
duction of C 2 . An example of reverse annealing 
in the experiments of T. Nakano, et al. is shown 
in Figure 16. These curves are not atypical. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
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